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Abstract 

We study how the influence of the shock wave appears in neutrino oscihations and the neutrino 
spectrum using density profile of adiabatic explosion model of a core-collapse supernova which is 
calculated in an implicit Lagrangian code for general relativistic spherical hydrodynamics. We 
calculate expected event rates of neutrino detection at Super-Kamiokande and SNO for various ^13 
values and both normal and inverted hierarchies. The predicted event rates of and Vf, depend on 
the mixing angle ^13 for the inverted and normal mass hierarchies, respectively, and the influence 
of the shock wave appears for about 2 - 8 s when sin'^ 2^13 is larger than 10"'^. These neutrino 
signals for the shock wave propagation is decreased by < 30% for Vf. in inverted hierarchy (SK) or 
by < 15% for Vf. in normal hierarchy (SNO) compared with the case without shock. The obtained 
ratio of the total event for high-energy neutrinos (20 MeV "^^^ MeV) to low-energy neutrinos 
(5 MeV < £^,y ^ 20 MeV) is consistent with the previous studies in schematic semi-analytic or other 
hydrodynamic models of the shock propagation. The time dependence of the calculated ratio of 
the event rates of high-energy neutrinos to the event rates of low-energy neutrinos is a very useful 
observable which is sensitive to ^13 and mass hierarchies. Namely, time-dependent ratio shows 
clearer signal of the shock wave propagation that exhibits remarkable decrease by at most factor 
~ 2 for in inverted hierarchy (SK), whereas it exhibits smaller change by ~ 10% for fe in normal 
hierarchy (SNO). Observing time-dependent high-energy to low-energy ratio of the neutrino events 
thus would provide a piece of very useful information to constrain ^13 and mass hierarchy, and 
eventually help understanding the propagation how the shock wave propagates inside the star. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

About twenty events of supernova neutrinos were detected from SN1987A . Some 
neutrino features such as neutrino temperatures and the energy carried out by neutrinos that 
lad been obtained from the data of SN1987A were consistent with theoretical expectation 
3|-|6|. However, the detected number of the neutrino events was still too small to find out 
more details about other flavors of the supernova neutrinos and the explosion mechanism. 
Although, we have not yet obtained neutrino events from the second core-collapse supernova 
which is the rare event of the century, detection of next supernova neutrinos is highly 
desirable to obtain important information on the neutrinos and the explosion mechanism. 
The expected information from the supernova neutrino observations is classified into two 
categories of implication for neutrino physics and for supernova physics. 

Large- volume underground detectors are now operating to detect various neutrino events. 
For example, Super-Kamiokande (SK), which is at the Kamioka mine in Japan, is a water 
Cherenkov detector, filled with 50,000 ton pure water (32,000 ton fiducial volume for the 

n 

burst mode and 22,500 ton for the other modes) [7[. KamLAND is a liquid scintillator 



detector with 1,000 ton active volume (8|. Sudbury Neutrino Observatory (SNO) has oper- 
ated as a heavy water Cherenkov detector, filled with 1,000 ton fiducial volume ^. Now, 



a new generation experiment, SNO+, is planned to be constructed [10|. We achieved re- 
markable development of the neutrino physics, especially about neutrino oscillation from 
the solar, atmospheric and reactor neutrino experiments with these detectors. In addition, 
the supernova neutrinos are thought to be fascinating targets of these detectors. If one 
supernova explodes at the Galactic center, thousands of neutrino events will be observed 

ll|. Furthermore, megaton-size neutrino detectors will detect xlO^ 



in Super-Kamiokande 



supernova neutrino events in the future 
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16(1 . The neutrinos are released directly from 



the core, in which the extreme physical condition of very high density is realized. Therefore, 
the supernova neutrinos would become the probe of such an environment. 

Neutrino flavor change by the neutrino oscillations relates to neutrino oscillation param- 
eters, i.e., the mixing angles, mass hierarchy, and CP violation phase. Most of neutrino 



oscillation parameters have been determined by the various neutrino experiments |17H19|. 
However, the mixing angle 613 is not precisely constrained, and only the upper bound is 
known from reactor experiments (e.g., 20|). In addition, the neutrino mass hierarchy, i.e.. 
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the sign of Am 



2 

13 



TTln 



ml, and the CP violation phase remain unknown. However, 



among many studies about imphcation on these unknown neutrino parameters from future 
supernova neutrinos 
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29| . there are several proposed possibilities that the detection of 



the neutrinos from the next Galactic supernova would constrain the neutrino oscillation 
parameters and identify the mass hierarchy more precisely 
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29|. 



Most of the supernova neutrinos are released for about 10 seconds after the core bounce, 
and interact with electrons when the neutrinos propagate through stellar matter. There- 
fore, Mikheyev-Smirnov- Wolf stein (MSW) effect is to be taken into account in the neutrino 



oscillations of the supernova neutrinos. The MSW effec 



; on the supernova neutrino signal 
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30 



32|). For example, the initial 



has been investigated previously in literatures (e.g., 
progenitor mass dependence of the early neutrino burst taking account of the MSW effect 
was studied in 



fied in 



34 



33| . and the Earth matter effects of the supernova neutrinos were identi- 



37j . Moreover, neutrino spin- flavor conversion in supernova has been studied 



considering numerically in 38l-l40| and analytically in 4ll, |42 1 . 

Recently, the effect of the shock wave on MSW effect of neutrino oscillations in supernova 



was studied 
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43| . This effect appears as a decrease of the average energy o 



Uf. in the 



case of normal mass hierarchy (or z/g in the case of inverted mass hierarchy) [4j, |45|]. Since 
the shock wave passes through H-resonance region, whose typical resonance density is ~ 
10^ g/cm'^ in a few seconds after core bounce, the adiabaticity of the resonance changes. 
The time dependence of the neutrino events monitors the density profile of the supernova 



26-28 
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46| . Therefore, using MSW effects embedded in the supernova neutrinos, we 



are able to find the density profile of the supernova. However, in many previous studies 
the influences of the shock wave on the neutrino events were discussed schematically using 
simplified and parameterized shock- wave profiles |27l |. 

There are many supernova simulations taking account of various explosion mechanisms. 
However, there are a few studies in which the neutrino oscillations were calculated by using 
simulation results because almost all supernova simulations were performed in only core 
region although MSW effect occurs at outer envelope of the star. In practice, it is still 
very hard to carry out multidimensional simulation with neutrino transport to proceed to 
the shock propagation from the core throughout the envelope. Therefore, it is important 
for the studies of the neutrino oscillations to calculate the shock propagation throughout a 
supernova not only in the core but also in the envelope. 
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In this paper, we study how the shock wave propagation and the neutrino oscillation pa- 
rameters affect the supernova neutrinos. We calculate neutrino oscillations of the supernova 
neutrinos by using supernova simulation result in adiabatic explosion model. In our model, 
we calculate all processes of core collapse, bounce, and shock propagation continuously in a 
unified manner. In order to clarify detailed dependence of the result on the neutrino oscilla- 
tion parameters and the shock wave, we calculate survival probabilities, energy spectra, and 
the event rates of three flavors of supernova neutrinos to be detected with SK and SNO. 

In Sec. [ni we introduce our numerical method. In Sec. Illlt we show time evolution of 
the survival probability and energy spectra of the supernova neutrinos. We also study the 
event rates. We analyze the dependence of these neutrino signals on neutrino oscillation 
parameters and the shock effect. We discuss the ratio of high- to low-energy neutrino events 
which will be detected by SK and SNO. We discuss the difference of the expected supernova 
neutrinos using both parameterized density and the density obtained by the simulations of 
supernova explosion in Sec. IIIIEI We also discuss other effects of supernova explosions and 
supernova neutrinos which should affect the time evolution of the neutrino signal. Finally, 
we conclude our paper in Sec. HVl 



II. NUMERICAL METHOD 



A. Supernova Model 

In order to obtain the detailed density profile in a supernova explosion and to examine 
the MSW effect on supernova neutrinos, we use one dimensional simulation result of the 
supernova. We model the supernova exp losion using an implicit Lagrangian code for general 
relativistic spherical hydrodynamics 47|. T he numerical tables of Shen's relativistic equation 
of state (EOS) 481] and Timmes's EOS [43] are adopted for the high and low density matters 
in th. code, respectively. We adopt the presupentova model of a 15M, sta. provided by 
Woosley and Weaver (WW95) [50] as the initial condition. We use the distribution of the 
electron fraction in WW95. We perform calculations of all processes throughout the 
core collapse, bounce, and shock propagation continuously by adiabatic collapse with fixed 
electron fraction to follow the shock wave for a long time scale (> 10 s) and a wide density 
range, consistently [51|]. We thus succeeded in the calculation of the propagation of the 
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shock, which is generated by the collapse of the iron core, until the shock wave reaches 
the hydrogen rich envelope. Throughout one calculation, one finds that the density of the 
central core reaches ~ 10^^ g/cm'^ and the shock wave propagation for beyond the hydrogen 



envelope (~ 1 g/cm^) is solved simultaneously 52l. |53|. 

Figure [T] shows the calculated density profiles as a function of radius for every 1 s from 
to 7 s after the core bounce. The shock wave reaches the oxygen-rich layer in 1 s and passes 
through the helium layer in 10 s. The horizontal lines show the density at the H-resonance 
of the neutrino energy. The shock wave reaches the H-resonance in about 2 s. The density 
behind the shock wave hardly decreases because we neglected neutrino cooling of the proto- 
neutron star in the present study. Therefore, the shock velocity is kept almost constant at 
9 X 10® cm/s in this model. When the neutrino cooling of the proto-neutron star is taken 
into consideration, the density and the pressure behind the shock would inevitably decrease. 
As a result, the shock wave would be decelerated. We study detailed analysis of the cooling 
effect of the proto-neutron star elsewhere. 

We assume the neutrino energy spectra with Fermi-Dirac distributions. Numerical simu- 
lations of supernova explosions which include neutrino transport suggested that the resultant 
neutrino spectra at low-to-intermediate energies do not exactly follow Fermi-Dirac distribu- 
tions. However, since we are interested in the event number of neutrinos to be detected in a 
water Cherenkov detector, the neutrino signals are sensitive to high energy neutrinos because 
the neutrino-induced reaction cross sections in the detector are proportional to the neutrino 
energy. As large as the neutrino energies above ~ 10 MeV are concerned, therefore, Fermi- 
Dirac distributions with finite chemical potentials may be justified as one of the reasonable 
approximations to the results of numerical simulations of the neutrino transport ( 5J] ; and 
references therein). Although several simulations suggest different neutrino temperatures 
from one another, they all satisfy a common hierarchy T^^ < Tp^ < Ty^. The interactions 
between v^. and neutrons freeze out at lower temperature and density than the interactions 
between z/g and protons do inside the neutralized supernova core, while the interactions be- 
tween fx (= t'/x, t'/x, '^T, and Pr) and nucleons freeze out at much higher temperature and 
density because interact with matter only through the neutral current. In the present 
study, the neutrino temperatures are set to be T^^ = 3.5 MeV, Tp^ = 4 MeV, and T^,^ = 7 
MeV associated with the neutrino chemical potentials fi^e = 7.28 MeV, fip^ = 10 MeV, and 
/ij,^ = MeV. These parameter values are well approximated to the time integrated spectra 
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of the neutrinos obtained from a supernova simulation in Livermore Group ll|]. Previous 
studies of energy spectra of supernova neutrinos in Galactic chemical evolution of ^^B and 
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56| and the r-process nucleosynthesis [57|, l58|] indicated that the neutrino temper- 



atures are constrained by the theoretical fit to the observed elemental abundances of these 
nuclei : T^^ = 2.5 - 3.5 MeV, Tp^ = 4 - 5 MeV, and T^^ =4-7 MeV, still satisfying 
lierarchy condition T^^ < Tp^ < T^^, where the uncertainty due to finite chemical potentials 



55| are taken into account in the inferred temperature values. Our adopted values T^^ = 3.5 



MeV, = 4 MeV, and T^^ = 7 MeV are in reasonable agreement with these constraints. 
The energy spectrum of each flavor of neutrinos is presumed to be independent of time. 
The total neutrino energy carried by neutrinos is set equal to 3 x 10^^ erg and the energy is 
equipartitioned to three-flavors of neutrinos and anti-neutrinos. We assume the exponential 
decay of the neutrino luminosity with a decay time of 3 s. 

B. Neutrino Oscillation 



In order to calculate neutrino oscillation effect on the neutrino spectrum, we solve the 
time evolution of the neutrino wave function along the density profile of our supernova 
model. The time evolution of the neutrino wave function is solved using the equation 



d 



dt 
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where AEij = Amjj/2Ei^, Amf^ = m| — m^, Gp, Ey and ne(r) are the energy difference 
between two mass eigenstates, mass squared differences, the Fermi constant, the neutrino 
energy, and the electron number density. In case of anti- neutrinos, the sign of \/2GF'ne 
changes. U is the Cabibbo-Kobayashi-Masukawa (CKM) matrix 



u 



C12C13 



■5l2Cl3 



-S12C23 - Ci2S23Si3e*'^ C12C23 - Si2S23Si3e"^ 



Size 
S23C13 



v 



— Cl2C23Sl3e*^ — C12S23 — Si2C23Si3e*^ C23C13 



(2) 



where Sij = sin 9ij, and Cij = cos 9ij (i 7^ j = 1,2,3) [59|. We here put the CP violating 
phase 6 equal to zero in the CKM matrix. 
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The neutrino oscillation parameters are taken from the analysis of various neutrino ex- 
periments 0, Q, except for 6*13 19 1, as 

sin^ 2^12 = 0.84, sin^ 2^23 = 1-00, (3) 

and 

Am?2 = 8.1 X 10"^eV^ jAm^g] = 2.2 x lO^Wl (4) 

We calculate the neutrino survival probabilities for both normal (Am^g > 0) and inverted 
(Am^3 < 0) mass hierarchies for four cases of sin^ 2^^i3 value: sin^ 2^13 = 10~^, 10^^, 10~^ 
and 10~^ 

The energy spectra of the neutrinos passed through the exploding supernova, 
are obtained from the spectra of the neutrino emitted from the neutrino sphere, (j)°^^{Ej,), 
multiplied by the survival probability in accordance with the following equation. 



/,SN 

kSN 



PiE,) 
P{E,) 
1 - PiE,) 1 - PiE,) 2 



1 - P{E,) ^ 
1 - P{E.) 
P{E,)+P{E,) J 



( 



iorg 



{E.\ 



\ 



mEu) 

\<l>l?iEu) J 



(5) 



+ (t>u.r + + (Ppt)^ P and P are survival probabilities of Ue and u^, 



where = ^{ 
respectively. 

Neutrinos change largely their flavors at resonances of neutrino oscillations. There are 
two resonances for supernova neutrinos (e.g., 2J]). The resonance points of higher and 
lower electron number densities correspond to H-resonance and L-resonance, respectively. 
The electron number density at the resonance point is written as 

1 Am^ 



Up 



cos 26*, 

E,, 



(6) 



2V2GF 

where Am^ and 6 correspond to [Amfgl and ^^13 at H-resonance and to Am^g and 612 at 
L-resonance, respectively. The flavor change at a resonance strongly depends on adiabaticity 
of the resonance. When the resonance is adiabatic, large flavor change occurs. The neutrino 
adiabaticity is estimated by 7, 

Am^ sin^ 29 n. 
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2E^ cos 2^ 1^1 



(7) 



If 7 is larger than 1, the resonance becomes adiabatic, but if it is small, resonance becomes 
non-adiabatic. The change of adiabaticity mainly occurs at H-resonance owing to the shock 
propagation. When the shock wave reaches the H-resonance, the density gradient becomes 
large and the value of 7 decreases. This leads to a change of survival probability of neutrinos. 
Since the resonance point is inversely proportional to the neutrino energy, this effect appears 
from low-energy side and moves toward high-energy side according to the shock wave prop- 
agation. As shock wave propagates outward, the density at the shock front decreases and 
the resonance condition is satisfied for higher energy neutrinos. 



C. Detection of Supernova Neutrinos 



The expected event number of the neutrinos in a water Cherenkov detector can be ex- 
pressed as 



dEjt - ^'^^ "^^^^^ M dKdt ""^^"^ M 



(8) 



where N is the detection number of neutrinos, is the energy of electron or positron, Ey 
is the energy of neutrino, Ntar is the target number, ri(^Ee) is the efficiency of the detector, d 



is the distance from the supernova. 



dE^dt 

section function of the neutrino energy 



is the neutrino spectrum, and (y{Eu) is the cross 



59|. We assumed a supernova at the center of 



Milky Way ((i=10 kpc), and we neglected the Earth effect 60|]. We assume the detection at 
SK and SNO. 

The mass of SK is 32000 ton and the solvent is pure water. SK mainly detects the Ve 
events with the reaction z/g + p — > e+ + n. The event number is obtained by integrating over 
the angular distribution of the events. As for the efficiency of the detector, we assumed as 
follows for simplicity: r]{Ee) = for E^. < 7 MeV and f](Ee) = 1 for Ef.> 7 MeV in accordance 



with SK Phase II 



6l| . The finite energy resolution of the detector was neglected here. The 



z/g and z/e events are evaluated taking account of the following four reactions. 



Ue+P 






+ n 


Ve + e~ 






+ e- 


z7e + e~ 






+ e~ 


+ e" 






+ e- 



(9) 
(10) 

(11) 
(12) 
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The energies of positron and electron are evaluated as Ee+ = Ey + nip — nin — for 
z7g + p — > n + e"*" and Ee = — ^ for + e~ — )■ z/ + e~, respectively. The cross sections of 
the above four reactions are adopted from jsQl]. 

The mass of SNO is 1000 ton and the solvent is heavy water. SNO detects not only 
events but also z/g events with charged and neutral current reactions. As for the efficiency of 
the detector, we assumed as follows for simplicity: ri(^Ee) = for E^ < 5 MeV and ri(^Ee) = 1 



for E,> 5 MeV 



62|. We evaluate the z/g and event numbers using the four reactions, 



z/g + c? — )■ p + p + e~ (13) 

z/g + rf n + n + e"^ (14) 

u + d p + n + u (15) 

z/ + e" z/ + e". (16) 



The main reactions of Ue and z/g are the first and second reactions, respectively [63|] . The 
cross sections of these reactions are adopted from |64]. The cross sections become larger as 
the neutrino energy increases. Although the SNO experiment has already ended, we also 
evaluate the neutrino events by SNO because the detection efficiency of z/g is high and similar 
in next generation detector like SNO. 



III. RESULT AND DISCUSSIONS 



A. Survival probability 

We calculate the survival probabilities of z^g and z/g from the supernova, P{Ey) and P{Ey). 
Figure [2] shows the calculated results. Left panels of Figure [2] are P{Ei,) in the case of normal 
hierarchy in every 1 s, and right panels of Figure [2] are P{Ey) in the case of inverted hierarchy. 
Red, green, blue, purple, sky blue and orange lines correspond to the results in t = 0, 1, 2, 3, 
4 and 5 s after core bounce. We showed the results for the values of sin^ 2^13 = 10^^, 10^^, 
10~^ and 10~^ from top to bottom. In the case of sin^ 2^13 = 10~^ and 10^'^, P{Ey) and 
P{Ey) are about when shock wave does not reach H-resonance (t = — 2 s). However, 
in later time (t = 2 — 5 s), P{Ey) and P{Ey) become finite. These effects appear from 
low-energy side and gradually shift toward high-energy side as the time passes by. In the 
case of sin^ 2^13 = 10~^, P{Ey) and P{Ey) hardly change with time, and are about 0.3 and 
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0.7, respectively. In the case of sin^ 26*13 = 10~^, P{Eu) and P{Eu) are about 0.2 and 0.4 
in t = — 2 s, respectively. These survival probabilities indicate intermediate adiabaticity 
between the cases for large ^13 and small 6*13. At the later times they become large as shown 
in the large 6*13 case. 

We can understand these behaviors considering the shock wave propagation. The adi- 
abaticity of H- resonance is estimated by 7 in Eq. ([7]). When the shock wave reaches 
H-resonance region, the adiabaticity changes. If sin^ 26*13 is large and there is not the shock 
wave at the resonance region, 7 is larger than 1, and the resonance is adiabatic. As a result, 
the survival probabilities are 0. When the shock wave reaches a resonance, large density 
gradient reduces 7 and the resonance becomes non-adiabatic. Thus, the survival probabili- 
ties become finite during the shock passage in the resonance region. In contrast, if sin^ 26*13 
is very small, 7 is smaller than 1, and the resonance is non-adiabatic. Consequently, the 
survival probabilities are finite, and the influence of the shock wave hardly appears to the 
survival probability, since H-resonance totally in non-adiabatic regardless the existence of 
the shock wave. In the case of sin^ 2^13 = 10~^, however, the influence of the shock wave is 
not clearly seen though the sin^ 26*13 is large. This is because that sin^ 26*13 is very large and 
7 is larger than 1 even when |^^| becomes large. 

7 can be rewritten by the ratio of the length characterizing the neutrino oscillations 
Lose = x^¥^, where Am^ = {{2^20 p E^n^{r) / Am'^ - cos2^)2 + sin^ 2^}i/2Am2, to the 
length of the level crossing region, i.e. the level crossing length 5r = sirm| j_*ie 1-1 

o o o ; 00 cos 20 1 We dr ' 

resonance. Therefore, the adiabaticity of the neutrinos can be evaluated by Lose and 5r 
59| . The resonance is adiabatic when the level crossing length is larger than the oscillation 
length. Lose ^ <^'"- Figure [3] shows the level crossing region and the oscillation length of 
neutrinos {Ey = 5 MeV) in the case of inverted mass hierarchy as a function of the radius. 
The left panels of Figure [3] are at t = s, and the right panels are at t = 3 s, respectively. 
Shown are the results for sin^ 2^13 = 10~^, 10~^, 10~^ and 10~^ from top to bottom. Red and 
green lines display the oscillation lengths of Ue and z/g, and blue line shows the level crossing 
length using ^13. The oscillation length of z/g becomes extremely large on the resonance, 
and its behavior depends on sin^ 2^13 below or above the resonance strongly. On the other 
hand, the oscillation length of z/g changes rather gently as a function of radius, and strong 
dependence on sin^ 2^13 is not clearly seen. As sin^ 26^13 becomes smaller or the shock wave 
reaches the resonance, the level crossing length becomes shorter. 
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In the left panels (t = s), the level crossing length at sin^ 2^13 = 10~^ and 10~^ are larger 
than the oscillation length at all region. Therefore, the resonance is adiabatic. On the other 
hand, in the case of sin^ 2^13 = 10~^ and 10~^, the level crossing length at the resonance is 
almost same as or slightly smaller than the oscillation length of z/g. Therefore, the resonance 
is non-adiabatic. In the right panels (t = 3 s) when the shock wave reaches the H-resonance, 
the level crossing length, 6r, is of the same order or smaller than the oscillation length Lose 
for all cases of the mixing angle, expect for sin^ 26'i3 = 10~^. This satisfies non-adiabatic 
condition 7 < 1 as discussed above. We understand, therefore, that the resonance becomes 
non-adiabatic by the effect of the shock wave. In Figure Ej the influence of the shock wave 
appears about 2 s after core bounce. This result is consistent with our simulation (see Figure 

CD. 

In the case of normal hierarchy, the survival probability of Ug does not change much 
because there is no resonance in the u sector. Therefore, the survival probability of in 
normal hierarchy is always ~ 0.7 regardless of the value of sin^ 2^13 or independently of the 
shock wave. 

We note that there is L-resonance in u sector even in the case of inverted hierarchy. 
However, the value of 612 which is related to L-resonance is very large (see Eq. ([3])), and 
the level crossing length of our simulation is not as small as the oscillation length at the 
resonance. Therefore, 7 at the L-resonance is larger than 1. As a result, the survival 
probability of z/g in inverted hierarchy does not change drastically, and stays always ~ 0.3. 

B. Supernova neutrino spectrum 

We calculate the supernova neutrino spectra using the survival probabilities. Figure H] 
shows the spectra of Ug. Left panels of Figure H] show the spectra in the case of normal 
hierarchy, and right panels show the results of inverted hierarchy. Figure [5] is same as Figure 
Hlbut for Uf, spectra. Red solid and blue dotted lines are the spectra with and without shock 
wave, respectively. In order to clearly observe the shock wave effects, we display these ratio, 
0with/0 without 5 of the spectra with to without shock in lower part of each panel. 

We see clearly the shock wave effects in the spectra in Figure H] in normal hierarchy. 
In the case of sin^2^i3 = 10"'^ and 10"'', an enhancement in low energy component of the 
neutrino spectra is seen when the shock wave reaches the H-resonance. At later times, the 
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influence of tlie sliock wave on tlie spectra moves from the low-energy side to the high-energy 
side. The effect at later times is seen as a reduction of the high energy component of the 
neutrino spectrum. The influence of the shock wave in the case of sin 26*13 = 10^2 (Joes not 
clearly appear in spectra for the same reason as discussed in the previous section. In the 
case of sin^ 26*13 = 10~^, the neutrino spectra with shock are not different from the spectra 
without shock. H-resonance is non-adiabatic even without the shock wave. We do not see 
any shock effects on the spectra in inverted hierarchy. 

From Figure |5l the effect of shock wave appears in the spectra of z/g in the case of inverted 
hierarchy. The dependence of the Ve spectra on sin^ 26*13 is almost the same as the dependence 
of the Vf, spectra in normal hierarchy. 

The spectrum of v^. at the surface of the supernova is written from Eq.dS]) as 



The initial spectra are 0°^^ > 0°^^ in low energy side, and (0°^^ — 0™^) is positive. Therefore, 
0°!^'^ increases when P is not 0. On the other hand, the initial spectra are 0™^ < 0°^^ in 
high energy side, and (0°^^ — 0°^^) is negative. Therefore, 0°|j"^ decreases when P is not 0. 
These increase and decrease appear in the observed neutrino spectra as the result from the 
influence of the shock wave. In the lower part of each panel, we can see this effect clearly. 
We note that the energy that satisfies 0°^^ = 0°^^ is about 18.7 MeV, and the energy that 
satisfies 0^'^ = 0°^s is about 23.3 MeV. 

C. Event rate of supernova neutrino 

We discuss here the predicted event rates of the supernova neutrinos to be detected with 
SK in the manner described in section III CI The upper parts of each panel of Figure [6] show 
the expected energy-integrated event rates of v^. 



where the integrand in right hand side is the expected event number of neutrinos defined 
by Eq. ([8]), and ii^th are the threshold energies which are equal to 7 MeV and 5 MeV for 
SK and SNO detectors, respectively, as explained in section III CI Left and right panels of 
Figure E] are for the cases of normal and inverted hierarchies, respectively. The lower part of 



0^f(E.) = P(E.){0-s(E.)-0-/(E.)} + 0-/(E.). 



(17) 




(18) 
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each panel shows the ratio of the event rates with and without shock. We calculate for the 
four cases of the mixing angles: sin^ 2^^i3 = 10^^ (red), 10~^ (green), 10~^ (blue) and 10~^ 
(purple). 

In the case of normal hierarchy (left panel of Figure [H]), the event rate of Vf. of each 
parameter is not separated clearly from one another. The event rate is commonly about 70 
at if: = s and decreases to about 3 at t = 10 s. On the other hand, the ratio of the event 
rates of z/g with and without shock shows some different 6'i3 dependence (left lower panel of 
Figure |6]). After 2 s, we can see slight enhancement of this ratio. We find from this change 
that the shock front does not reach H-resonance before ~ 2 s. Once the shock wave reaches 
H-resonance region the event rate increases from the low-energy neutrinos because of energy 
dependence of the resonance density Eq. (|H]). This ratio then decreases after 3 s because the 
event rate of the high-energy neutrinos decreases when the shock wave propagates through 
the resonance. The obtained time profile thus depends on the ^13 values. The effect of 
shock propagation appears most clearly in the case of sin^ 2^^i3 = 10^'^ (green). However, 
the decrease in the ratio of the event rates is by at most 15 % in the case of small 6'i3. It is 
smaller than the case for sin^ 2^13 = 10^'^ at any time for the other 61^ values. 

In the case of the inverted mass hierarchy (right panel of Figure |6]), the time variation 
of the event rate of z/g does not depend on 6'i3. Shock wave effect is not clearly seen (lower 
part of right panel of Figure |6]). The adiabaticity of z/g is not influenced by the shock wave 
because z/g is not related to H-resonance in the case of inverted hierarchy. Therefore, the 
spectra of do not change for any values of 61^. 

Figure [7] is the same as Figure El but for z/g. In the case of the normal hierarchy (left panel 
of Figure [7]), the event rate and the shock wave effect of z/g are insensitive to 6^13 throughout 
the time profile, and event rate decreases from about 3,200 at t = s to about 100 at t = 10 
s. This is attributed to the fact that the adiabaticity of Ve is not influenced by or the 
shock wave as discussed in the previous section. 

In the case of the inverted hierarchy (right panel of Figure [7j), the behavior of event rate 
of z^g is different from z/g in normal hierarchy. At t = s, the event rate is about 4,200 for 
sin2 2ei3 = 10-2 and lO'^, and it is about 3,200-3,500 for sin2 2^i3 = 10"^ - 10"^ The 
event rate for large ^13 (i.e. sin^ 2^13 = 10"^ and 10~^) is larger than that for small 6'i3 (i.e. 
sin2 2^i3 = 10~^ and 10~^). When H-resonance is adiabatic, almost all z/g are completely 
converted from z/^.. On the other hand, when H-resonance is non-adiabatic, some z/g remain 
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as z/g even after the neutrinos get out of the star. As the result, the fraction of v^. converted 
from z/a; is small. The average energies of the three flavor neutrinos just emitted from the 
proto-neutron star satisfy the following relation; E^^ < Ep^ < E^,^. Therefore, the average 
energy of z/g in adiabatic case becomes higher than that in non-adiabatic case. As the main 
reaction cross section for is proportional to E"^ as discussed in section III C\ the detection 
probability of high-energy neutrinos is larger than that of low-energy neutrinos. Therefore, 
the event rate of z/g for large ^13 (adiabatic case) becomes much larger than that for small 
^13 (non-adiabatic case). This is the reason for large split between the expected event rates 
for large ^13 (i.e. sin^ 2^13 = lO^^ and 10"^) and small ^13 (i.e. sin^ 2^13 = 10"^ and 10"^). 
We find difference of event rates by ~ 1, 000 at t = between the two cases. 

Furthermore, the ratio of the event rates for z/g which exhibits the shock wave propagation 
effect changes dramatically when sin^ 2^13 is large (right lower panel of Figure [7]). The ratio 
of z/g increases at around 2 s, and decreases after 3 s, which is similar to z/g event rates in 
normal hierarchy (left lower panel of Figure [6]). The mechanism is quite similar to each 
other so that the shock wave propagation through the H-resonance starts from low energy 
neutrinos to high energy neutrinos, which eventually decreases the event rates. However, the 
decrease in the ratio of the z>e event rates is remarkable by at most 30 % for sin^ 2^13 = 10~^. 
On the other hand, the ratio changes weakly for the other 613 values because H-resonance 
is non-adiabatic regardless of the shock wave, although its degree is still larger than the z/g 
event rate (see left lower panel of Figure [6]). 

The event rate of Ue is extremely larger than that of z/g (for a given neutrino oscillation 
parameter set). This is mainly because the total cross section of z^g-induced reactions is about 
10^ times larger than that of z/g. Moreover, the total cross section of z>g-induced reactions 
is almost proportional to the square of neutrino energy, but the total cross section of z/g- 
induced reactions is linearly proportional to the neutrino energy. Owing to this different 
energy dependence of the cross sections, one can expect that the effect of the shock wave 
should be identified more clearly in the z/g events. 

In practice, since the shock effect is folded in the observed event rates of z/g and z/g, 
we should carry out careful theoretical analysis. The observed event rates are statistically 
uncertain of the order of their square root. As for z/g in normal hierarchy, the statistical 
error is the same order as the change in the ratio ~ 15% of the event rates with and without 
the shock effect because absolute number of expected event rates are small ^ ^ 70 (Figure 
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[6]): Actually, ^ ~ 10 at t ~ 6 s when the shock effect becomes maximum strength, and 
thereby \J^§-/^ ~ 0.3. On the other hand, for in inverted hierarchy, the statistical error 
is expected to be small enough because of large absolute number of expected event rates 
^ < 4000 (Figure CD: In this case, since ^ ~ 500 at t ~ 6 s, sj^/^ ~ 0.05 which is 
much smaller than the maximal change in the ratio ~ 30% of the event rates. 



D. Ratio of high- to low-energy neutrino events 

There is a potential difficulty in identifying the influence of the shock wave upon the 
observed neutrino events because time evolution of the original neutrino spectra, which 
is unaffected by the shock wave as well as the MSW effect, is unknown. The original 
neutrino spectra must be assumed theoretically a priori. We therefore look for another 
useful observable that should show a clear signature of the shock wave effect even when we 
do not know the time evolution of the original neutrino spectra. 



1. Time-integrated ratio of neutrino events 

We consider now the ratio of the high-energy component to the low-energy component 



of time- integrated neutrino events. This ratio -Rx is defined by [32|, 

»10s /-eOMeV ^2jY 

dEAt 



dEfdt ^ 



T3 — JOs J20McV "-^e"-'' q\ 

— ^los ^20MeV ^2jY ' ^^^^ 

-dEpdt 



Os 



where x refers to SK or SNO and -Eth = 7MeV (SK) and 5 MeV (SNO) as explained below 
Eq. ( |T8l) . We set here the boundary between high- and low-energy components to be 20 
MeV because the original energy spectra satisfy the conditions 0°^^ = and 0p^^ = 
at E^ = 18.7 MeV and 23.3 MeV, respectively. 

Figure [8] shows Rsk vs -Rsno for various mixing angles sin^ 26'i3 = 10~^(red), 10"'^ (green), 
10~^(blue) and 10~^(purple). Left panel of Figure [H] is in normal hierarchy, and right panel 
is in inverted hierarchy. Closed and open circles are the calculated results with and without 
shock wave effect, respectively. Note that the scale of each panel is different. 

In the inverted hierarchy, we see large variations of the ratios; -Rsk ~ 1.8 - 4.6, and 
-Rsno ~ 2.3 - 4.3. In addition, Rsk and -Rsno show a clear correlation because the most 
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dominated reactions for SK and SNO detectors are the z/g-induced charged current reactions; 
Ve+V e'^ + n for SK and v^ + d n + n + e^ ioi SNO. Both ratios Rsk and -Rsno increases 
with the mixing angle ^13. For a given ^13, the ratios in the case with shock are smaller 
than those without the shock. On the other hand, in the normal hierarchy, the ratio shows 
a small dependence on sin^ 2^13 for i?sK and -Rsno- In practice, we see only small variation 
of the ratios; -Rsk ~ 1-76 - 1.8, and -Rsno ~ 2.2 - 2.5. These characteristics in the parameter 



32|. 



dependence of -Rsk and -Rsno are consistent with the result of ref. 

The parameter dependence which we confirmed in -Rsk and -Rsno as discussed above could 
be an important observational signature to discriminate neutrino mass hierarchy, normal or 
inverted, as well as the mixing angle, ^13. In normal hierarchy, the ratios -Rsk and -Rsno are 
small regardless of the mixing angle 61^ and the shock effect. Furthermore, -Rsno changes 
a little according to the value of 6'i3 and -Rsk hardly changes. We defined -Rx as the total 
events of all flavors. Most of the events detected in SK is Vf.- As a result, the variation 
of Ve events is hindered by larger events unchanged by the variation. On the other 
hand, the events of SNO contain the events of Ve as much as that of Ve-, because the cross 
section of i/g-induced reaction is of the same order as that of z/g-induced reaction. Therefore, 
we can see the variation of v^.. However, the variation of v^. is small because the number 
of events is not so different between Ve and Ve- In inverted hierarchy, the ratios -Rsk and 
-Rsno correlate with the mixing angle 6*13 and the shock effect. Larger 6'i3 value indicates 
larger -Rsk and -Rsno ratios. However, the shock effect reduces -Rsk and -Rsno even in large 
^^13 value. Therefore, we do not distinguish the value of 6*13 and the shock effect in -Rsk and 
-Rsno- 



2. Time- dependent ratio of neutrino events 

In the previous subsection, we discussed the ratio of high-energy to low-energy events, 
-Rx, which is the ratio of the neutrino events integrated over the time. We also evaluated 
in Sec. IlllCl the time-dependent event rates integrated over the neutrino energy. Figures [6] 
and [7] indicate that some effects of shock wave propagation could be observed at later times 
when the shock wave reaches H-resonance. However, we apparently lose some important 
information on the effects of shock wave in either case. In this subsection, we explore for 
the signature of the shock wave effects by taking account of double differential as to both 
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time and energy. 

We here define time-dependent ratio of tlie events of high-energy to low-energy neutrinos, 



where x refers to SK or SNO. Figure [H] shows this ratio -Rx(^) as a function of time. The 
values of sin^2^i3 are 10~^, 10"'^, 10~^ and 10^^ from top to bottom panels. Left panels 
are for normal hierarchy and right panels are for inverted hierarchy. Solid and dashed lines 
are the calculated results with and without shock wave, respectively. Note that the scale of 
each panel is different. 

First, we show -Rsno(^) in normal hierarchy in left panels of Figure [9l The effect of shock 
wave is seen, but it is a very small effect of at most about 10 %. As for observed v^. spectrum 
in normal hierarchy, 70 % is the original and 30 % is the original Vf. in the non-adiabatic 
case, though almost 100 % is the original v^^^ in the adiabatic case. On the other hand, 70 % of 
the observed z/g is the original v^, and 30 % is the original in inverted hierarchy. Moreover, 
the shock effect of Vf. in normal hierarchy is washed out by neutral current and charged 
current of Ve because the target of neutrino detection is deuteron. Therefore, the variation 
of the event rate ratio is small in the normal hierarchy. There is, however, possibility of 
finding shock wave effect if much larger events are detected by the next generation detector. 

Second, we show RsK{t) in inverted hierarchy in right panels of Figure [91 The event rate 
ratio RsK{t) is constant in early times and strongly depends on sin^ 2^13. It is 4.7 and 1.8 in 
the case of sin^ 26*13 = 10^^ and 10~^, respectively, at t = s. The value of -RsK(i^) is large 
when sin^ 26*13 is large, because the H-resonance is adiabatic and the average energy of z/g is 
high. The small -Rsk(^) value is due to non-adiabatic state of H-resonance. 

When the shock wave effects are included and sin^ 2^13 is 10~^ (large), i?sK(^) greatly 
changes with time. Although -Rsk(^) is 4.7 until 2 s, it decreases to ~ 2, which is closer to 
the ratio in the case of non-adiabatic state in 4 - 8 s. This decrease is due to the change 
of the adiabaticity of H-resonance according to the shock propagation. The H-resonance 
is adiabatic before the shock arrival to the resonance. In contrast, the H-resonance is 
non-adiabatic while the shock wave propagates in the resonance, and RsK{t) approaches 
-Rsk(^) ~ 1-8 which is the ratio in small sin^ 26*13 case (see the bottom panel). 




(20) 
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When sin^ 2^13 is as small as 10~^, we see only small differences of RsK{t) between the 
cases with and without the shock wave and hardly see the time dependence. The H-resonance 
is non-adiabatic whenever the shock wave is on the resonance or not in this case of small 
sin^ 2^13. 

To conclude, we can see a clear signature of the shock wave in the time-dependent high- 
energy to low-energy ratio of -Rx(^) although the influence of the shock wave is less clearly 
seen in the ratio of time integrated event rates, R^- Therefore, the observations of the time 
evolution of -Rx(^) is quite an important observable to constrain the neutrino oscillation 
parameters and to find the effects of shock propagation in supernovae. 

E. Further discussion 

Fogli et al.js^] found qualitatively similar result to ours in the calculated event rates 
with and without forward shock in the case of inverted hierarchy for sin^ 2^13 = 10~^. Their 
results are however different from ours quantitatively for several reasons. For one, their event 
rates at the neutrino energy 45 ± 5 MeV are calculated by assuming the next generation 
detector of 0.4 Mton pure water that is 12.5 bigger than 32 kton for SK. We can scale our 
calculated results by 12.5 times in order to remove this apparent difference. The second 
reason is that they adopted completely different supernova model from ours. Although it is 
hard to reconcile the difference between the two models, we tried to compare at least the 
effect of shock wave propagation in the following manner. 

We should, first, refer to our result of the event numbers at the same neutrino energy 45 
± 5 MeV at 5 s after the core bounce for sin^ 2^13 = 10^'^. This aims to compare the two 
results of the event rates so that the shock wave effect appears most remarkably in both 
calculations in the case of inverted hierarchy. Second, we should remove the bias based on 
different supernova models especially arising from different density profiles from each other. 
For this purpose we normalize our event number without shock to Fogli's event number 
without shock. We thus obtain finally the following event numbers; 135 with shock and 
300 without shock in inverted hierarchy. These scaled and normalized numbers should be 
compared with Fogli's results; 130 with shock and 300 without shock. Note again that 
the latter number in either calculation is the same by definition of normalization. Our 
numerical calculations show 55.0 % decrease in the event number as a resultant net effect of 
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shock wave propagation, which is in reasonable agreement with 56.7 % decrease as shown 
by their calculations. We thus conclude that their theoretical calculations and ours of the 
neutrino event rates agree with each other despite several essential differences: We solved 
in Eq. (1) both flavor conversion of neutrinos and supernova density profile numerically in 
our method described in Sec. 1, while Fogli et al. j^/i applied an artificial model of density 




profile of shock wave propagation without solving dynamical supernova explosion. 

For further comparison, let us analytically estimate the decrease in the neutrino events 
in the non-adiabatic case (with shock) compared with the adiabatic case (without shock). 
In the non-adiabatic case, we define the survival probability of i>e as P = 0.7 in our model, 
and we use the result of z/g in the case of normal hierarchy in Fogli's model. We set P = 
in the adiabatic case. Thus estimated analytical results show 62.3 % and 63.0 % decrease 
in the neutrino events from adiabatic to non-adiabatic cases in our model and their model, 
respectively. These values are close enough to each other, and ~ 60% are not very different 
from ~ 55.0% which was inferred from numerical calculations as discussed in the previous 
paragraph. This fact justifies that our treatment of analytical estimates are reasonable. 
Therefore, regardless of all possible differences between the two models, the supernova neu- 
trinos provide a powerful tool to indicate the shock wave propagation inside supernova if 
sin^ 26'i3 is large in the inverted hierarchy. 

We assumed that CP violating phase is zero and sin^ 2^23 = 1 throughout this study. 
However, CP-phase is unknown and there is an uncertainty in the mixing angle 623. Here 
we discuss the dependence of the shock effect on neutrino flavor transitions on CP-phase 
and 623- Transition probabilities of neutrinos in arbitrary density profile have been studied 
theoretically in [65|], where it was found that the transition probabilities of z/g — )■ z/g and 
Ue — )■ z^e do not depend on CP-phase and mixing angle 623 at all. Furthermore, one can 
show from their formula in [65:| that the sum of the transition probabilities of — )■ z/g and 
Ur — z/g (z/^ — )■ z/g and Ur — t'e) are totally free from CP-phase and 623. Since the spectra of 
Ufj, and Ur {v^ and v-r) are the same at the neutrino sphere, the transition probability from z/^ 
or z/t- (z/^ or 1)^) to z/g (z/g) can be written as a half of the sum of the transition probabilities 
of z/^ — )■ z/g and Vr ^ {i^fi and z/,- — )■ z/g). The transition probability from z/^ or z/,- 

(z/^ or z/t-) to z/g (z/g) does not depend on CP-phase and 6*23. We thus conclude that there is 
no effect of CP-phase and mixing angle 6'23 on the supernova neutrino spectra. 

If the influence of the shock wave is seen very early {t < Is) in the observation of the 




20 



supernova neutrinos, it might mean that the shock wave reaches H-resonance (~ 1000 g/cm'') 
very early after core bounce. Then, we would expect that this kind of supernova explosion 
exhibits a strong shock wave effect in the direction of observer. Moreover, assume that we 
could know viewing angle from the axis of a supernova in optical and radio observations. 
If the viewing angle is small and the influence of the shock wave is seen very early, this 
supernova has strong shock like a narrow beaming jet along the axis of supernova explosion. 
Therefore, we would expect that this supernova explodes in the mechanism associated with 
the jet formation 66|-|72| . If the viewing angle is large and the influence of the shock is seen 
very early, we would expect that this supernova has a rather wide jet and may be close to 
the spherical symmetric explosion. 

On the other hand, if the influence of the shock wave is seen in the supernova neutrinos 
at relatively later time, and the viewing angle is small, we would expect that this explosion 
is also close to the spherical symmetry without jet. Therefore, we would expect that this 
supernova explodes in the mechanism without a jet. Details of the influence of neutrino 



oscillation on jet explosion are discussed in [73[. 

Recently, oxygen emission-line profiles from supernovae are observed, which could be a 
signature of an aspheric explosion [tJ]. We would study the asphericity of the supernova 
explosions more in detail from such an optical observation and simultaneous detection of 
the supernova neutrinos. 

The flavor-exchange effects by the neutrino self-interactions might also be important 
because of their huge flux immediately after the emission out of the proto-neutron star (e.g., 
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76|). This could change initial neutrino spectrum from what we assumed here. 



IV. CONCLUSION 



We investigated the neutrino signal to study the effect of the shock wave propagation as 
well as the dependence on unknown neutrino oscillation parameters of mass hierarchy and 
^13 by carrying out numerical calculations of the neutrino event number using the simulation 
result of a supernova explosion. We followed adiabatic collapse of iron core, core bounce, and 
the shock wave propagation in the stellar envelope for a long time ( ^ 10 s) using general 
relativistic hydrodynamical code and realistic density profile. Then, we could calculate 
detailed time evolution of the event number rate and spectra of the supernova neutrinos. 



21 



The expected event rate of Pg in the case of the inverted mass hierarchy and that of Ve 
in the normal mass hierarchy depend on the mixing angle ^13. When sin^ 2^13 is larger than 
~ 10~^, the influence of the shock wave appears in the observation after 2 s in our model. 
It is the time when the shock wave reaches the H-rcsonance. Therefore, the shock effect and 
the constraint on 9i^, can be inferred even by the event rate integrated by the whole energy 
range. However, it is difficult to distinguish the influence of the shock wave and the neutrino 
oscillation parameters only from the time evolution of the event rate. 

The time-integrated ratio of the events of high-energy to low-energy neutrinos is another 
indicator to examine the shock effect and constrain mass hierarchy and ^13. The parame- 
ter dependence of the time integrated ratio -Rx (x = SK or SNO) could be an important 
observational signature to discriminate neutrino mass hierarchy, normal or inverted, as well 
as the mixing angle, ^13. Both of the ratios i?sK and i^sNO correlate with the value of ^13 
in inverted mass hierarchy. Therefore, the neutrino oscillation parameters ^13 and the mass 
hierarchy, might be able to be constrained by this ratio. The shock effect reduces the ratios 
i?SK and -RsNO- This effect is not distinguished with smallness of ^13. 

The time-dependent ratio of high- to low-energy neutrino events is more preferable indi- 
cator to find out the shock effect clearly. The ratio decreases after 2 s for -Rsno(^) in normal 
hierarchy and for RsK{t) in inverted hierarchy. The decrease is the most remarkable in the 
case of sin^ 2^13 = 10~^ and more clearly in inverted hierarchy. The dramatic decrease in 
the ratio could be a clear signal for the shock wave effect and would constrain the minimum 
value of sin^ 2^13. Therefore, observations of the time-dependent ratio of the high- to low- 
energy neutrino event are an important observable to find the effects of shock propagation 
in supernovae. 
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FIG. 1: The density profiles as a function of radius for every 1 second from to 7 seconds after 
the core bounce. The horizontal lines show the density at the H-resonance point of each neutrino 
energy. 
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FIG. 2: Calculated results of the survival probabilities, P[Ey) and P{Ey), in every 1 second. Left 
panels are the survival probabilities of v^, P{Eu), in normal hierarchy, and right panels are the 
survival probabilities of I'e, PiEu), in inverted hierarchy. Red, green, blue, purple, sky blue and 
orange lines are t = 0, 1, 2, 3, 4 and 5 sec after core bounce. The value of sin^ 2^13 is 10~^, 10~^, 
10~^ and 10^^ from top to bottom. 
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FIG. 3: The oscillation length, L, 



^ 2g ^j^^ i^^^i crossing length, 6r - iJ^^^^i-i 
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cos 2d ' Ue dr ' ' 

of neutrinos for 5 MeV in the case of inverted mass hierarchy as a function of a radius r in units 
of solar radius Rq = 6.96 x 10^'^ cm. The left panels are at t = sec, and the right panels are at t 
= 3 sec, respectively. We calculate by the value of sin^ 20i3 = 10~^, 10^'^, 10~^ and 10~^ from top 
to bottom. Red and green lines are the oscillation length for 13 mass eigenstate for neutrinos and 



anti-neutrinos, and blue line are the level crossing length. 
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FIG. 4: Upper part of each panels are Ve spectra ^{Vf^ in every 1 second. Solid (red) and dotted 
(blue) lines are the calculated results with and without shock wave, respectively. Lower part of 
each panel is the ratio of the spectra with shock to without shock, <\>with{Ev) j 4>w%thout{Ev)^ at 
t = l(red), 2(green), 3(blue), 4(purple) and 5(sky blue) sec. The value of sin^ 2^13 is 10~^, 10^'^, 
10~^ and 10^^ from top to bottom. Left and right panels are in normal and inverted hierarchy, 
respectively. 3q 
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FIG. 5: The same as Figure HJ but for spectra. 
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FIG. 6: Upper part of each panels shows the expected energy-integrated event rates of V(. with SK 
as a function of time, and lower part shows the ratio of these event rates with shock to those without 
shock for various mixing angles sin^ 2^13 = 10~^(red), lO^^(green), 10~^(glue) and 10~^(purple). 
Left and right panels are in normal and inverted hierarchy, respectively. 
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FIG. 7: The same as Figure [6l but for 
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FIG. 8: Calculated correlation between Rsk and -RsnOj which are the ratio of the time-integrated 
events of high-energy neutrinos (20 MeV < < 60 MeV) to those of low-energy neutrinos (5 
MeV < -Ejy < 20 MeV) for the SK and SNO detectors, respectively, for various mixing angles of 
sin^ 2^13 = 10^^ (red), 10^^ (green), 10~^(blue) and 10~^ (purple). Left panel is for normal hierarchy, 
and right panel for inverted hierarchy. Closed and open circles are the calculated results with and 
without shock wave, respectively, as denoted. Note that the scale of each panel is different, and a 
closed square in the right panel corresponds to a range in the left panel. 
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FIG. 9: Time-dependent ratios of the events of high-energy neutrinos (20 MeV < < 60 MeV) to 
those of low-energy neutrinos (5 MeV < < 20 MeV) for the SK and SNO detectors for various 
mixing angles of sin^ 2^13 = 10~^, 10~^, 10~^ and 10~^ from top to bottom panel. Left panels are 
the ratio iisNo(^) for normal hierarchy, and right panels are the ratio RsKit) in inverted hierarchy. 
Solid and dashed lines are the calculated results with and without shock wave, respectively. Note 
that the scale of each panel is different. 
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